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In this paper we analyze a recent experiment conducted in an anechoic chamber, where the 
scattering of microwaves from an array of metallic cylinders was measured. This is a system which 
' displays chaotic scattering in the short wave limit. The analysis of the experimental data is aimed 

' at elucidating the effects of the underlying chaos. We describe a robust numerical method which 

provides the scattering matrix for any number of non overlapping reflecting cylinders. We use 
this method to calculate the scattering from one, two and three cylinders, and to compare the 
results in the numerical simulations with those obtained in the experiment. Both simulations and 
measurements validate the presented theory. 
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I. INTRODUCTION 



> 

' Scattering of electromagnetic waves in the microwave region is governed to a large extent by the underlying ge- 
ometrical optics or ray dynamics, which provides the dynamical skeleton over which the residual interference and 
diffraction effects can be superimposed. In many instances, the ray dynamics is very complex, and it displays chaotic 
features which earns this phenomenon the name "chaotic scattering" . What makes the dynamics chaotic is that the 
scattering trajectories which come from infinity, are trapped for a long time within the interaction region before they 
escape to infinity. In the interaction domain, the motion is exponentially unstable, and therefore the outgoing trajec- 
tories hardly "remembers" the initial conditions. This results in very complex reaction patterns, fractal distributions, 
self similarity and other attributes which are typical of classical chaotic dynamics. As is well known, the complex 
S ' patterns in dynamical systems are due to the existence of a "strange attractor" in phase space. In chaotic scattering, 
: which is not a dissipative process, the "strange attractor" is replaced by a "strange repeller" which consists of all the 
(~| trajectories which are trapped forever in the interaction domain. The complexity of the ray (or trajectory) dynamics 
O manifests itself also in the corresponding wave phenomena (be it electromagnetic or acoustic waves or all possible 
■ de-Broglie -quantum- waves describing electrons, atoms or nuclei). The purpose of the present paper is to illustrate 
J> , the importance of "chaotic scattering" in the context of electromagnetic wave physics. 

The system which we studied in this work consists of three parallel metallic cylinders which scatter electromagnetic 
waves in the microwave region. This study has been performed by means of simulations and, what is more important, 
cd ' through measurements in the anechoic chamber of the European Microwave Signature Laboratory (EMSL) [Q . This 
system illustrates in an elegant and simple way the "generic" features of chaotic scattering. It is also a paradigm in 
the theoretical studies of this subject (see e.g., @ 0, |^ and references cited therein). 

In the experiment, we measured the scattered fields in the TM mode (with the electric field parallel to the cylinders 
axis). The experimental set-up was planned such that one could neglect the coupling to the other polarization 
component. The antennas were in the far-field region of the cylinders (at a distance of about 10 m) and therefore the 
wavefront can be considered to be planar. The incidence is normal to the cylinders axis. This reduces the problem to 
the treatment of a two-dimensional (2-D) scalar field scattering of three cylinders in a planar geometry. 

This paper is organized as follows. In Sect. |l[ we introduce the formulation of the problem and describe a robust 
numerical method which provides the scattering matrix for any number of non overlapping metallic cylinders. We 



use this method to calculate the scattering from one, two and three cylinders. Then, in Sect. [II, we review the 
most important features which characterize the chaotic ray dynamics in the present system. The identification of 
the fingerprints of chaotic scattering of the rays in the actual experiments, where the scattering of waves of finite 
wavelengths are measured, is addressed in Sect. For this purpose we make use of the Eikonel (semiclassical) 
approximation and propose a few correlation functions which are shown to depend explicitly on functions which 
characterize the chaotic dynamics of rays. In Sect. ^ we turn to the description of the measurement set-up and the 
analysis of the experimental data. We compare the results in the numerical simulations with those obtained in the 
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experiment. We show the actual correlation functions and compare them to the predictions of the short wave theory 
proposed. We end this work with a summary of our conclusions and a few suggestions for further experiments. 



II. FORMULATION OF THE PROBLEM 

Consider N parallel infinitely long conducting circular cylinders as shown in Fig. |^. The axes of the cylinders 
are parallel to the z-axis of a cylindrical coordinate system. The center of the jth cylinder of radius aj is located at 
{Rj, (pj). A monochromatic plane wave with the electric field parallel to the axes of the cylinders (TM mode) impinges 
at an angle 9inc- We assume a time dependence of the type exp{—iwt) . The scattered field is scalar and must satisfy 
the planar wave equation 

(V2 + fc2)* = (1) 

where k denotes the frequency wavenumber. This equation is solved subject to Dirichlet boundary condition on the 
boundaries of the cylinders (i.e. the total tangential electric field vanishes on the surface of all cylinders). Thus 

^-Wlr^r, =0 (2) 

where rj is a point on the boundary of the cylinder j. At large observation distances, the solution may be decomposed 
into a summation of incoming and outgoing cylindrical waves as follows, 

" (2^ ^ 
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which is used as a definition of the scattering matrix S. This matrix contains all the information needed for a complete 
description of the scattering process and can be calculated exactly by making use of the KKR method @. The various 
matrix elements Si^i' describe the scattering between states of different angular momenta I and Because of fiux 
conservation S must be unitary. 

Following the work of Gaspard and Rice |^ we proceed to the calculation of the S matrix for the general case of 
TV non-overlapping cylinders. To this end, one uses the Green's theorem in order to convert Eqs. and (^) into an 
integral equation of the form 

/ d'r' [vI/fe,(r')(Vr", +k2)Go(r,r') -Go(r,r')(V?, +k2)VI>k,,(r)] = (4) 

JD 

(f dr'h{v') [^fc,,(r')Vr'Go(r, r') - Go(r, r')Vr'*k,i(r')] • 

JdD 

Here, D is the domain exterior to the cylinders but inside a circle with a boundary dDoo which is large enough to 
enclose all the cylinders. dD is the boundary of this domain, and the normal pointing outside is denoted by fT-(r'). 
The boundary dD of D is the union of dooD and the boundaries of the A^— cylinders and djD, (j = 1, . . . , A^). The 
free (outgoing) Green function is given by 

Go(r,r') = -^<)(k|r-r'|) (5) 
where h'^^ is the zeroth-order Hankel function of the first kind. 



Choosing r on the j cylinder the first part of Eq. <^ vanishes and we get: 

N 



= Y.f ^^'"(^J') [G(rj,rj')Vr'*k,i(rj') 



(6) 



where the contour integral has A^-|- 1 parts. By calculating this integral we obtain the gradient of the wavefunction 
at the cylinder boundaries as 

A = CM-^ (7) 
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where the matrix A defines the gradient of the wave function 



n(rj)V*fc,i(rj) = Aij„,e™«J, 1, m = 0, ±1, ±2, ±3, ... (8) 

m— — oo 

and 9j is the angular coordinate of the point rj . 

The matrix M describes the multiple scattering between the A'^ cylinders and is found to be 

M,„,,„, = Sr^m'S.r + (1 - S,,-)^^^^^^^H^^l^,ikR,,,) Q,rim,m') (9) 

where Rj,j' — Rj' .j is the distance between the j and j' cylinder, and aj is the radius of the jth cylinder. Note that 
M has the structure of a KKR matrix and is the generalization of the results of Gaspard and Rice ^ . The matrix C, 
is equal to 

m ) = exp [i{m (j>ji_j — m'{(f)jji — tt))] , (10) 
0,,'(m,m')-(-ir""'C;',,(m',m) 

and contains-besides a phase factor-the angle aj/j of the ray from the center of the cylinder j to the center of the 
cylinder j' as measured in the local coordinate system of cylinder j. In the simple case of one cylinder M = 7r/2i. 
Finally, the matrix C is given by 



2» ^i^^ Ji-AkR' ) 



^1]-^ = —^ n-j , ■ (11) 



with {R^ , $j) denoting the polar coordinates of the center of the jth cylinder as measured in the global coordinate 
system. 

If we now take r outside the smallest circle which encloses all the cylinders, but inside the domain D, the first 
part of Eq. (^ provides 5'/c_i(r). We then obtain the S matrix at large distance from the scatterer by propagating 
the wave function, whose gradient on the cylinders is now known. At the last stage of the calculation we use the 
asymptotic expansions of the Bessel functions at large distances r and the S matrix is identified by comparing with 
the asymptotic expression (|^). The final result reads, 

5 = / - iCM-^D (12) 

with the matrix D defined: 

A,™ = -7ra,J™_j(fci?^Vj(fca,)e-™*'. (13) 

The method presented above has the advantage that the calculation of the S matrix is essentially reduced to the 
algebraic problem of inverting the matrix M defined in Eq. (^). In principle, this matrix is infinite. However, in 
order to carry on with numerical calculations, we need to approximate it by a finite square matrix of dimension L. 
The integer L is chosen in such a way that the truncated matrix is unitary to within the desired accuracy of the 
calculation ^. 

The normalized scattered field is given in terms of the scattering matrix S as 

/(0i,„0,,,;fc) = -|-— ^ ^ e-''^(<^-^-(^/^\Si,,-6,,,)e^^'('^-^-(-/^\ (14) 

^ ' 1,1' — — OD 

where 0inc and 9sct denote the incidence and observation azimuth angles, respectively. 
The 2-D bistatic radar cross section (RCS) is then simply computed as 

CTbist (6*1110, 6*801 ;fc) = \f{Ohic,Osct;k)\'^ ■ (15) 

The presented formulation takes multiple scattering to any order and is valid for any working frequency, provided 
that the number of cylindrical modes used in the scattering matrix L is sufficiently large. In the next section, we will 
focus on the chaotic behavior of the scattered fields in high frequency limit. 
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III. RAY DYNAMICS 



Under the approximation of the geometrical optics, the electromagnetic scattering is governed by the simple law of 
specular reflection at the boundary of the cylinders. Given an arbitrary frame of reference, any incoming ray can be 
identified by its direction 9i and by the nearest distance by which it approaches the origin - the impact parameter bi. 
The impact parameter is proportional to the angular momentum li = kbi where k is the wavenumber. In the domain 
of all possible initial rays {9i,bi), there is a compact set of values for which the corresponding rays will impinge on 
one of the cylinders. An initial ray in this relevant set will undergo multiple reflections by the cylinders, until it is 
scattered away emerging as the outgoing ray (0/, 6/). The function Of{bi; 9i) where 9i is kept fixed and bi varies within 
the relevant set is called the deflection function. 

For scattering on a single cylinder whose center is at the origin, we have 



The cylindrical symmetry implies bf — bi, and the actual deflection Of — 6i is independent of 9i. Adding another 
cylinder has two important consequences- the cylindrical symmetry is broken, and multiple scattering becomes possible 
in the vicinity of the trapped orbit which runs along the line which connects the two centers of the cylinders. Note 
that a ray which starts at infinity, can approach the trapped ray arbitrarily closely, but it must finally emerge. 
The resulting defiection function cannot be calculated analytically. Its complex structure is shown in Fig. |^a. The 
complexity is due to the fact that the trapped orbit is unstable. (The instability emerges because the reflections are 
induced by concave mirrors). The longer the scattering orbits dwells next to it, the more rapid are the fluctuations 
in the deflection function. The concept of "dwell time" is essential for the following discussion and it will be properly 
defined for a general scattering ray (independently of the number of cylinders it may encounter). Consider a scattering 
ray coming from 9i,bi, and after Q reflections it goes out with Of,bf. Denote the reflection points on the various 
cylinders by r^ , with 1 < j < N . Then the dwell time (measured in units of length) is defined as 



Here, c is the speed of light, fii and nj are unit vectors in the incoming and the outgoing directions. The last two terms 
above, remove the arbitrariness in the definition of the point at which the dwell-time stop-watch starts (and stops) 
to tick. Fig. ^ shows the dwell time which correspond to the deflection function for two cylinders. The correlation 
of the complexity of the deflection with the dwell time is evident. 

To generate chaotic scattering, one needs at least one more cylinder, which should be placed such that no cylinder 
interrupts (shadows) the lines which connect the boundaries of its neighbors. In this case, there exists a Cantor 
set of trapped orbits which can be encoded by the following symbolic dynamics: Denote the three cylinders by the 
letters A,B,C. Then, to any infinite string of letters {Xi}°^^, Xi G {A,B,C}, with Xi ^ X^+i there corresponds 
an unstable trapped orbit. If the string is N periodic {Xi — X^+at) for all i, the trapped unstable orbit is periodic. 
This is the "strange repeller" which is responsible for the chaotic dynamics in our system: A ray which comes from 
infinity and is scattered by the three cylinders is affected by an infinite set of trapped orbits. It can dwell next to 
any of them for some time, and then, due to the intrinsic instability, it can either be trapped next to another member 
of the repeller, or scatter out. The resulting dynamics is much more complex than what we saw in the two cylinders 
system, and it has the following features: 

• The deflection function is singular on a Cantor set of bi values, showing self similar structures which occur on 
all scales (see Fig. ||). 

• The dwell time function is singular on a Cantor set of bi values, which is correlated with the singular, self similar 
structures observed in the deflection function (see Fig. ||). 

• The deflection function displays also smooth sections, which are due to single scattering from the convex hull of 
the boundary. These trajectories correspond to short dwell times, and we shall refer to this component of the 
scattering system as "direct" . 

• Excluding the direct component, the rest of the scattering dynamics is ergodic - any bundle of neighbor incident 
trajectories emerge with (0/,6/) values which arc uniformly distributed over the relevant set defined above. 



6f{bi;e.i) = ei + 2a.rcsinbi for |bi| < 1. 



(16) 




(17) 
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The distribution of dwell times is asymptotically exponential (see Fig. 

-Pdwon(T) w 7rexp(-r7^) 



(18) 



and 

-fr = {l-dH)X (19) 

where dn is the Hausdorff dimension of the Cantor set of singular points on the bi axis, and A is the Lyapunov 
exponent which characterizes the instability of the strange repeller. 

• The distribution of impact parameter transfers (defined as the difference between the final and the initial impact 
parameters) for the ergodic component is approximately uniform 

/^impact lOj ~ I Q otherwise ' 

D is an effective diameter of the part of the scatterer which induces the ergodic scattering component (see 
Fig.|b). 

Much more can be written on the mathematics and physics of chaotic ray scattering. The interested reader can 
find more details and references in the articles cited above. The material presented above suffices for the purpose of 
the following discussion. 



IV. CORRELATION FUNCTIONS - THE FINGERPRINTS OF CHAOTIC DYNAMICS 

The scattering amplitude /(^i, 02', k) is the transition amplitude to scatter from 6i to 62 at a wavenumber k. In the 
Eikonel approximation the scattering amplitude can be expressed as a superposition of amplitudes. Each amplitude 
corresponds to a geometrical ray which is incoming from 61 and is out-going in the direction 62- These rays are 
identified in the following way. Consider the deflection function Of{hi \ Oi) where 9i is fixed at the value of interest 6*1 
and one looks for the impact parameters 5| which satisfy the equation 

e2 = ef{ht-e,). (21) 

The geometrical rays which are incoming at (6|, 9i) emerge at 62, and hence they correspond to the desired transition 
6*1 -^62- Due to the fact that the deflection function oscillates wildly as a function of hi, there is an infinite set of 
values hf which satisfies (^). The Eikonel approximation for the scattering amplitude is 

/(0i,02;fc)«(27rfc)-^^ 

s 

Here, the sum goes over all the rays which satisfy (|2l|), and the dwell time for each ray ( [L7| ) appears explicitly, and 
it is considered as a function of the parameters 9i , 62 which define the transition. The pre-exponential factors can 
be interpreted as the partial scattering amplitudes brought by the individual rays. As a matter of fact, the sum of 
their absolute squares is the cross section in the geometrical optics approximation. It is obtained by deleting from 
the Eikonel cross section all the terms which come from interferences of contributing amplitudes. The integers Us are 
the Maslov indices. 

The scattering amplitude gets contributions not only from the multiply reflected rays, but also from the "direct" 
component of rays which scatter from the convex hull of the cylinders. The statistical treatment we propose is not 
applicable for this component, and we have to subtract it from the measured or calculated data. This is done using 
the following reasoning. The "direct" component corresponds to very short dwell times, and therefore its dependence 
on k is very smooth. Therefore, if one subtract the k smoothed scattering amplitude one remains with the desired 
statistical component. The smoothing is done over a k interval of the size Afc centered about k. In the present 
case, the smoothing interval was taken to be the entire interval where the measurement was carried out. Thus, the 
"statistical" component of the scattering amplitude is 



deid92 



(22) 
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k+Ak/2 



/(0i,02;fc)dfc 



(23) 



The f{Oi, 02', k) is the primary object of our statistical study. Its Fourier transform with respect to k provides the 
spectrum of dwell-times for the rays which support the transition — > (?2- When averaged over the incoming and 
outgoing directions, it should coincide with the expression (nq) i.e. 



/ 



dk f {91,02; k)e- 



(24) 



A comparison between the experimental measurement and the classically expected exponential decay is shown in 
Fig. ||. The behavior is clearly exponential over some 7 decades. The slope of the curve is equal to jr ~ 1.068 ns~^, 
in reasonable agreement with the value 7,- = 0.95 ns^^, calculated from the ray dynamics. 
We will now define the angle averaged correlation function 



da 

— ( /(a, a + AO; k)r{a + e, a + e + A0; k) 

ZTT 



(25) 



Here, we consider a situation where the outgoing angle relative to the incoming direction A9 is kept constant, and one 
correlates the scattering amplitudes for different orientations of the target. The k averaging is the same as defined in 
(p3|). Using the Eikonel approximation (^) and the relations 



dT{ei,e2) 



- -bi 

= 62 



we can approximate the correlation function by 



C.(e) 



(26) 



(27) 



To derive this expression, we assumed that e is small so that the phase differences can be approximated by the 
leading order in their Taylor expansion. Off-diagonal contributions are discarded due to the k smoothing, and the 
pre-exponential factors are expressed in terms of the geometrical-optics partial cross sections 



Ps{a,A0;k) = 2iik 



861862 



with 61 ^ a 



a + A6. 



(28) 



The sum is extended over all rays s. We can now collect together all the partial cross sections which are due to rays 
for which the impact parameter transfer [llf^ — h')^'^) equals b 



(29) 



With this definition we get 



Ce(e) w ( / dbP,mpact{b; k)e \ 
\j —00 I k 



(30) 



Since we subtracted away the "direct" contribution, we can use the ergodicity of the chaotic component and estimate 
P{b) by the expression (^ ). 

In a similar way we can define the cross correlation 



Xe{e)^ I ^(^f{a,a + Ai6;k)r{a + e,a + e + A26;k)'^^ 



(31) 
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where, in contrast with (??) we correlate scattering amphtudes measured at two different detector configurations. 
The rays which support the scattering in the two configurations are completely unrelated, and therefore we expect 

Xg{e) = 0. (32) 

Another useful correlation function can be defined by 

Ct{i^)^ ^{f{a,a + M-k + ^)~r{a,a + M-k~^^))^ . (33) 

By repeating a similar argument as above, one can express Ck{£,) as a Fourier transform of the dwell time distribution 

/oc 
drPd^e«(r)e-^'^" . (34) 
-oo 

The dwell time distribution Pdweii{T) was discussed above, and the expression ( p^ implies that the |Cfc(i^)|^ is a 
Lorentzian with a width which is proportional to 77-. 

The two correlation functions obtained above are typical for scattering systems for which the resonances are overlap- 
ping and which display therefore Erickson type fluctuations . One of the most important achievements of the theory 
of chaotic wave scattering was to recognize the fact that Erickson fluctuations are the hallmark of the underlying 
chaotic ray dynamics. 

Several other correlation functions (e.g., Hanbury-Brown Twiss type correlations) can also be computed jj]. How- 
ever, the measurement requires the use of two receiving antennas at close (angular) proximity. This experiment was 
not yet carried out. 



V. MEASUREMENT DESCRIPTION AND DATA ANALYSIS 



Before comparing our numerical simulations with the experimental measurements we shall give a short description 
of the experimental setup. The measurements were performed in the anechoic chamber of the European Microwave 
Signature Laboratory. This chamber has a hemispherical shape with a diameter of about 20 m. It is equipped with 
two separate TX/RX antenna modules, which can be moved independently along a circular arch in a vertical plane, 
as well as 27 receiving antennas which are on the uniformly distributed hemispherical dome. The range from the 
center of the chamber to all the antennas is the same: 9.56 m. The system is fully polarimetric and operates in the 
stepped frequency mode. 

The measurement set-up used in the experimental validation is shown is Fig. |7| The target (consisting of one, two 
or three identical copper cylinders of height 30 cm and diameter 7.9 cm) was placed on a rotating table whose axis of 
rotation is parallel to the axes of the cylinders. The center of gravity of the three-cylinders coincides with the rotation 
axis (vertical axis). Two different antenna configurations were used, and they are referred to as Tl and T2, with 

• Tl configuration: TX at angle 89.5°; RX at angle 25.85° 

• T2 configuration: TX at angle -89.5°; RX at angle 25.85°. 

For each configuration, the bistatic scattered fields complex amplitude was measured at 36 equally spaced azimuth 
positions of the table at 10° intervals. The data were sampled within the frequency interval 4 — 20 GHz with a 
frequency step of 11.25 MHz. The acquired data in the frequency domain were empty room subtracted and gated in 
the time domain in order to isolate the response of the cylinders from the residual antennas coupling and eventual 
spurious reflections in the chamber. Then a bistatic calibration using a metallic sphere of diameter 30 cm placed at 
the focal point of the chamber, was applied in the HH and VV polarizations. 

Some examples of experimental and numerical radar cross section vs. the frequency v are shown in Figs. ^ and ^ 
for antenna configurations Tl and T2 and for one, two and three cylinders. As we see, the agreement of our numerical 
calculations with the experimental measurements is quite good within the entire frequency range. 

We invested some effort to understand the deviations which appear for the Tl configuration at the higher frequency 
range. In principle, the scattering of a single cylinder should not depend on the table orientation. However, the 
experimental data show distinct periodic modulation with peak value of w 15 %. This can only be explained as due 
to the finite geometry: the effective scattering angle does depend on the orientation since the cylinder was not placed 
at the center of the rotating table. 
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Finally, Fig. |To| shows the variations of the bistatic cross section for some specific values of frequency, when the 
azimuth position of the rotating table is varied. Again we are able to see a quite good agreement with the experimental 
results which lead us to the conclusion that our numerical calculations simulate quite well the experimental data. 

We turn now to the statistical analysis of the scattering matrix. More specifically we calculate the various correlation 
functions defined in section |^ in order to test the statistical predictions. Fig. |ll| displays the functions Ce^e) (??) 
and Xg{€) (??). The slight deviation of the cross correlation from the expected value is due to the fact that we 
perform a finite averaging over the frequencies. The autocorrelation function, however is consistent with Eq. (|30|). 

Another interesting quantity we can evaluate is the energy correlation function Cfc(x) defined in Eq. (??). From 
our semiclassical consideration, we expect it to have a Lorentzian form with half- width approximately equal to 7^. 
Our results, reported in Fig. are in agreement with the semiclassical prediction. 

A direct consequence of the above result, and of the underlying ergodicity is that not only the amplitudes but also 
the cross sections show similar correlations in energy. We expect 

where o'bist(fc) = o'bist(fc) ~ (o'bist(fc))fc is the fluctuating part of the cross section. Again, we calculate Ca{i') and 
recovered a satisfactory agreement with the semiclassical predictions. The results are shown in Fig. pT 



VI. CONCLUSIONS 



We have presented an algorithm that takes multiple scattering to any order and we can use it to simulate any 
experimental configuration. We present results which show that chaotic scattering can be unambiguously identified 
in the experimental data. Moreover, the various theoretical simplifications which were assumed, as well as the finite 
angular resolution in the experiment, did not blur the signatures of the chaotic ray dynamics. 

The main conclusion from the present calculations is that in the range of frequencies used in the experiments, 
the effects due to chaotic scattering are prominent in the present system and, in general in most systems in which 
electromagnetic radiation is scattered from metallic objects. The scattered radiation field is so complex, that it shows 
many features which would be typically attributed to a random field. 

One may investigate the subject further by performing e.g., the measurements of the Hanburry-Brown Twiss effect, 
or by considering other scattering systems such as cavities (see 0) or waveguide networks (see [^). We also believe 
that the observations made here can be used, and should be taken into account when heterodyn electromagnetic 
radiation experiments analyzed. 
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FIG. 2. (a) The deflection function for the two cyhnder system as a function of the imj^t parameter; (b) The corresponding 
dwell time as a function of the impact parameter 
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FIG. 3. The deflection function for the three cylinder system as a function of the imlpq,ct parameter. Three successive 
magnifications are shown, (the magnified domain is framed by dashes) and the fractal nature of the deflection function is 
obvious. 
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FIG. 4. (a) The deflection function for the three cyUnder for the entire range of impact pabpieter; (b) The correspond 
dwell time exhibits singularities which are correlated with the self-similar structure observed in the deflection function. 
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FIG. 5. (a) The probability distribution of the dwell times. The slope of the straight line is O.^^J^ ^CWl^^^ distribution 
of impact parameters for the three cylinder system. 
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FIG. 7. Sketch of the measurement set-up used in the experimental validation. 
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FIG. 8. Experimental (circles) and numerical (stars) bistatic cross section vs. fr§j^u*iVg[|ftJr2;|ie Tl configuration. 
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FIG. 9. Experimental (circles) and numerical (stars) bistatic cross section vs. B(;fec(J^p|c(7T^ltlie T2 configuration. 
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Figure 6 
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